Weak compressible magnetohydrodynamic turbulence in the solar corona 
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This Letter presents a calculation of the power spectra of weakly turbulent Alfven waves and fast magne- 
tosonic waves ("fast waves") in low-(3 plasmas. It is shown that three-wave interactions transfer energy to high- 
frequency fast waves and, to a lesser extent, high-frequency Alfven waves. MHD turbulence is thus a promising 
mechanism for producing the high-frequency waves needed to explain the anisotropic heating of minor ions in 
the solar corona. 
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The heating of the solar corona is a long-standing prob- 
lem. Measurements taken with the Ultraviolet Coronagraph 
Spectrometer (UVCS) have provided important constraints on 
coronal heating, showing, for example, that T± 3> 7]| and 
7j_ ~ 10 8 K for +5 ions at a heliocentric distance r of 
roughly two solar radii, where T± and 7j| are the temperatures 
for random motions perpendicular and parallel to the mag- 
netic field B. 1 1, 2] These measurements imply that the aver- 
age magnetic moment kgT^/B of +5 ions increases rapidly 
with r and show that +5 ions are heated by electric and 
magnetic-field fluctuations with frequencies CO comparable to 
or greater than the ions' cyclotron frequency £2. (If CO <C £2, 
the average magnetic moment is almost exactly conserved.) 

Different sources have been proposed for the required high- 
frequencywaves, including reconnection events in the coro- 
nal base heat-flux-driven plasma instabilities |7], 
and magnetohydrodynamic (MHD) turbulence |& |<J 1 1 C11 - 
An apparent difficulty with this last source is the finding 
that in incompressible and weakly compressible MHD turbu- 
lence there is little or no cascade of energy to high frequen- 
cies. However, incompressible and weakly 
compressible MHD neglect the fast magnetosonic wave ("fast 
wave"). In this Letter, a weak turbulence calculation is used 
to show that when fast waves are accounted for, MHD turbu- 
lence in low-p plasmas transfers energy to high-frequency fast 
waves and, to a lesser extent, high-frequency Alfven waves. 
(In the corona, p = %%p/B 2 ~ 0.01, where p is the pressure.) 
The high-frequency waves produced by MHD turbulence may 
be of importance not only for coronal heating, but for particle 
acceleration in solar flares as well. 11611 

The MHD momentum and induction equations with Lapla- 
cian viscosity and resistivity are 

fdv _\ _/ B 2 \ BVB - „ 

and 

= Vx (vxB)+r|V 2 B, (2) 

at 

where p is the density, v is the velocity, p is the pressure, 
and B is the magnetic field. In this Letter, the magnetic field 
is taken to consist of a uniform background field and a small- 
amplitude fluctuating field: B = Bqz + 5S. The pressure is 



discarded since p is taken to be <C 1 . The spatial Fourier trans- 
forms of v and b = 5B/y/4np can be written 

Vk = v a ,k e u ,k + v/,k k± + v St kZ (3) 

and 

bk = b a ,k e a ,k + b/,k e f . k , (4) 

where e a j = z x k± is the Alfven-wave polarization vector at 
wave vector k, k± — k±/k±, k± — k — k z %, and ef± — e a ^ x 
k/k. The Alfven-wave amplitude is given by = v a k^b a ^, 
and (since p <C 1) the fast-wave amplitude is given by ft 1 = 
Vf^±bf£- Upon neglecting nonlinear terms in the momen- 
tum and induction equations, one obtains dat/dt = ±.ik z v&aj; 
and dfc /dt = ±ikvAf k ± , where va = fio/V^p is the Alfven 
speed. Alfven waves have frequency ^fk z VA and propagate 
along magnetic field lines. Fast waves have frequency ^fkvA 
and can propagate in any direction. The v s ^z term in equa- 
tion (0 corresponds to the slow magnetosonic wave, which 
has a frequency that approaches zero as p — * 0. 

Weak turbulence consists of waves whose amplitudes are 
sufficiently small that nonlinear interactions between waves 
can be treated as a small perturbation to a wave's linear be- 
havior. Weak turbulence theory is based on the assumptions 
of random wave phases and approximately Gaussian statis- 
tics. 11711 These assumptions are problematic for acoustic tur- 
bulence, because sound waves propagating non-dispersively 
in the same direction interact coherently for long times. 1 1 71 
The same issue arises for fast waves, because of their 
acoustic-like dispersion relation. However, fast-wave interac- 
tions with Alfven waves and slow magnetosonic waves limit 
the interaction time for pure fast-wave interactions, which 
may allow weak turbulence theory to apply to MHD at low P 
even if it does not apply to acoustic turbulence. Although this 
issue remains unresolved, weak turbulence theory is a valu- 
able starting point for this difficult problem. 

To simplify the analysis, the slow magnetosonic wave is 
neglected, the density p is taken to be a constant, and (to 
maintain energy conservation when p is held constant) the 
v • Vv term in equation ([0 is replaced with \A ■ Vv, where 
is the part of the velocity associated with Alfven waves. A 
different approach was taken by El, who included slow 
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waves but neglected three-wave interactions that did not in- 
volve slow waves. Further work including all the nonlineari- 
ties is needed. The Alfven-wave and fast-wave power spectra 
for homogeneous turbulence are defined through the equations 
(4(4)*) =A±6(*-*i), and (jf (/*)*> = F±6(*- *i), 
where (...) denotes an ensemble average. It is assumed that 
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A+ = Al = A k , that F+ = Fj7 = F k , and that <a±/±) = 
( a kfk\) = ^- Rotational symmetry about the z axis is also 
assumed, so thatA^ = A(kx,k Z) t) and F k = F(kj_,k z ,t). Tak- 
ing the small-v and small-r| limits and employing the stan- 
dard weak-turbulence approximations, one obtains the wave 
kinetic equations, 



dt 



dt 



8va 



8va 



d 3 pd 3 qh{k p q) {8(q z )8( P± nl) 2 A q (A p ~A k ) + 8{k z + Pz + q)M l [M 2 F q {A p -A k )+M 3 A p (F q -A k )] 
8(k z + p z - q)M 4 [M 5 F q {A p -A k ) + M 3 A p {F q -A k )} + 8(k z +p-q) M 6 [M 7 F q {F p -A k )+M % F p {F q -A k )]} 



(5) 



d 3 p d 3 q8(k-p-q){9 sin 2 B[8(k-p- q)kqF p (F q -F k ) + 8{k + p- q) (k 2 F p F q + kpF q F k - kqF p F k )} 



+ 8{k-p z + q z )M 9 [M 10 A q (A p -F k )+M n A p (A q -F k )] + 8{k- p z - q)M 12 [M n F q (A p - F k ) +M l4 A p (F q ~ F k )] 
+ 8{k+p z -q)M l5 [M l6 F q (A p -F k )+M 17 A p {F q -F k )]}, 
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where 

M 2 = — p±m — (cosa+ l/2)(k±l + p±m + q±n), 

Mi = 2k±l + 2p±m + q±n, 

M5 = — pi_m+ (cosa— 1 /2)(kj_l + pj_m + qj_n), 

M-i = k±J(cosa- l/2)+p_Jn/2 +sinan(2p— q/2), 

Ma = fcjJ(cosy+ 1/2) - sin\)fm(2q- p/2) — q_Ji/2, 

Mio = p±m+ (cos0+ l/2)(kj_l + pj_m + qj_n), 

Mn = qi_n+ (— cos0+ 1 /2)(kj_l + pj_m + qj_n), 

Mn = sin0l(— k + 2q) +p±m(cosQ — cosa) + sinan(<7 — 2k), 

Mu = smQl(k/2-2q)+p A m(-cosQ+ l/2)-q A n/2, 

M16 = sin0/(— k + 2q) +p±m(cosa — cos9) + sinan(<7 — 2k), 

Mn = sinQl(k/2 -2q)+ p A m(cosB+ l/2)-q 1 n/2, 

Mi = M2 +M3, M4 = Ms +M3, Me = M 7 +Mg, Mg = Mm + 
Mn, M12 =Mi3 +M14, and M 15 =M 16 +Mi 7 . The quanti- 
ties 0, \|/, and a are the angles between z and the wave vec- 
tors k, p, and q, respectively. In the triangle with sides of 
lengths k±, p±, and q±, the interior angles opposite the sides 
of length k±, p±, and q± are denoted <3 k , G p , and <3 q , and 
/ = coso^, m = cosGp, n = cosa q , I = sinc^, m = sino p , 
and n = sina ci - The above form of the wave kinetic equa- 
tion makes use of the identities k± cos(a q — G p ) = q±n + p±m 
and k± sin(c q — O p ) — q±n — p±m. 

The right-hand sides of equations (0 and (the "colli- 
sion integrals") represent the effects of resonant three-wave 
interactions. The integrals sum over all possible wavenumber 
triads, while the delta functions restrict the sum to triads sat- 
isfying the resonance conditions k = p + q and <x> k = (x) p + (£> q , 
where co^ is the frequency at wavenumber k. The equa- 
tions Mi = M2 +M3, M4 = M 5 +M3, etc imply that dA k /dt 
(or dF k /dt) is positive at any wave number at which A k (or 
F k ) vanishes, provided the spectra are positive at other wave 
numbers. The wave kinetic equations thus ensure that the 
spectra remain positive (realizability). Since the dissipative 



terms have not been included, equations l|5} and (Jfjji conserve 
the energy per unit mass E — J d 3 k(A k +F k )/2 and have an 
equipartition solution F k =A k = constant. 

The 8(q z ) in the collision integral of equation l|5} is equiva- 
lent to 2vA$(k z VA — P z Va +<1zVa) and represents the frequency- 
matching condition for resonant interactions involving three 
Alfven waves ("AAA interactions"). The part of the colli- 
sion integral that contains this 8(q z ) is the same as the colli- 
sion integral for AAA interactions in incompressible MHD. 
This term represents interactions between oppositely directed 
Alfven waves, in which the field-line displacements caused 
by Alfven waves travelling in one direction along the mag- 
netic field [represented by A(q±,q z = 0)] distort Alfven wave 
packets travelling in the opposite direction. At k z = 0, only the 
AAA terms contribute to the right-hand side of equation 10, 
and the steady-state solution A(k±,k z = 0) « kj 3 can be ob- 
tained analytically with the use of a Zakharov transformation, 
as in the incompressible case 0. When A(k±,k z = 0) °c k^ 3 , 
and when non-AAA interactions are neglected, a Zakharov 
transformation yields A k k^ 3 for any k z . It can be seen from 
equation (|3J that the time scale Xa for AAA interactions to 
transfer Alfven-wave energy from k± to 2k± is determined 
by A(k±,k z = 0) and is independent of k z , consistent with 
physical descriptions of the Alfven-wave cascade, omul 
If the Alfven-wave energy per unit mass (8v rms ) 2 is dom- 
inated by wavenumbers of order some characteristic wave 
number ko, if the spectrum is quasi-isotropic at k ~ ko, and 
if A(k±,k z = 0) oe k^ 3 for k± > ko, then x A ^ v A]\k \ (8vrm S ) 2 l 
for &x ^> as in the incompressible case. 

The terms in the collision integral of equation © propor- 
tional to 8{k — p — q) and 8(k + p — q) represent resonant 
three-wave interactions involving only fast waves ("FFF in- 
teractions"). As can be seen from the delta functions, FFF 
interactions occur only when k is parallel or anti-parallel to 
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both p and q, indicating that these interactions transfer energy 
radially in £-space. The FFF terms are the same as the colli- 
sion integral for weak acoustic turbulence flill . up to an over- 
all multiplicative factor proportional to sin 2 0, and represent 
a weak form of wave steepening. As sin0 — * 0, the acoustic- 
like FFF interactions become less efficient because the fast 
waves become less compressive. If the non-FFF terms are ne- 
glected, then a Zakharov transformation can be used to show 
that Fk = cig(8)fc -7 / 2 is a steady-state solution to equation l|6) 
for any function g(0). When F k = cig(0)AT 7 / 2 , the energy 
flux in FFF interactions per unit mass per unit solid angle in 
£-space, e, can be obtained in the same way as for weak acous- 
tic turbulence fiUl . and is given by e = 9% 2 c\ sin 2 Bg 2 c2/16vA, 
where c 2 = £ 'dxln(l +x)[x(l +x)]- 5 / 2 [(l + x) 9 l 2 -x 9 l 2 - 
1] ~ 26.2. If the fast waves were forced isotropically and non- 
FFF interactions were ignored, then g — 1/ sin 8 in steady state 
so that e is independent of 0. The time scale if for FFF in- 
teractions to transfer fast-wave energy from k to 2k can be 
estimated by dividing the fast-wave energy per unit solid an- 
gle between k and 2k by the energy flux e. Ignoring numer- 
ical coefficients, one obtains if ~ va/[ci sin 2 9g(0)A; 1 / 2 ] for 
F k = cig(Q)k-y 2 . If the fast -wave energy were dominated by 
wavenumbers of order some characteristic wave number ko, 
with Fk = ag(Q)k~ 7 / 2 for k > ko, then c\ ~ (Sv^^ ) 2 k X J 2 , 
where (Svmis.f ) 2 is the energy per unit mass in fast waves. In 
this case, i F ~ v A /[(8v nnS: /7) 2 (A;oA;) 1 / 2 sin 2 eg(0)] for k > k Q . 

The terms in equations and l|6) containing Mi 
through M\i correspond to three-wave interactions involv- 
ing both Alfven waves and fast waves ("AAF and AFF in- 
teractions"). Such interactions exchange energy between fast 
waves and Alfven waves within resonant wavenumber triads. 
At small 0, the frequencies of fast waves and Alfven waves 
are comparable, and AAF and AFF interactions are efficient. 
For example, if F% ~ cifc _7 / 2 sin _1 and A k <C F k at small 0, 
then when 9 <gc 1 the largest contribution to dF^/dt comes 
from the term proportional to M\ ^F q F k and is —F k /iAf, where 
^af = (15va sin0)/(237t 2 cife 1//2 ) to lowest order in 0, a time 
scale that is <C if . The energy lost by fast waves in this case is 
transferred primarily to Alfven waves at the same wavenum- 
ber through the term in equation containing M%F p F q . If 
A k grows until A k = F# at small 0, then the term containing 
MijF q Fk is cancelled by the term containing M\^A p Fk to low- 
est order in 0, largely stemming the loss of fast-wave energy. 
AAF and AFF interactions thus act to make A k ~F k at small 0. 
However, the constant-energy-flux solution A k ~ F k oc k-y 2 is 
unsustainable, because as k increases energy is lost from the 
small-0 part of £-space to high k± through AAA interactions 
faster than it is replenished from small k by FFF interactions 
(if °c k~ x l 2 , 1a « k~ l ). The energy flux in FFF interactions 
at small must thus decrease with k as fast-wave energy is 
drained into Alfven waves and then transferred out to large k±. 
This process causes Fk to steepen relative to AT 7 / 2 at small 0, 
and results in Alfven-wave energy at \k z \ ^S> ko. On the other 
hand, for > 45°, the frequencies of Alfven waves and fast 
waves differ considerably, and AAF and AFF interactions are 



unable to make ~ F k at k 3> ko. In this part of £-space, AAA 
and FFF interactions dominate the right-hand sides of equa- 
tions (0 and (|6}, so that Fk « k~ n l 2 and A k = /z(£ z )&f[ 3 within 
the inertial range, where h(k z ) is some (decreasing) function 
of\k z \. 

To obtain quantitative solutions for Ak and Fk, equations l|5) 
and are integrated forward in time numerically with ini- 
tial spectra A k = F k = fc 2 exp(— k 2 /kfy. The isotropic forc- 
ing term C3& 2 exp(— k 2 jkfy is added to the right-hand sides of 
both equations (0 and (|6). The dissipation terms —ct,k 2 A k and 
—CA,k 2 F k are added to the right-hand sides of equations (0 and 
(|6) respectively, with 04 chosen so that in steady state dissipa- 
tion truncates the spectra at a wavenumber that is 3> ko- The 
numerical method conserves energy to machine accuracy in 
the absence of dissipation and forcing and will be described 
in a future publication. Steady-state spectra at late times are 
plotted in figures ^ an d E] and are consistent with the quali- 
tative picture described above. The Alfven-wave spectra are 
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FIG. 1: Alfven-wave power spectrum as a function of k±_ at differ- 
ent k z . 

°= h[k^jk~^ for k± \k z \ within the inertial range. Although 
h(k z ) decreases with increasing \k z \, there is some Alfven- 
wave energy at \k z \ 3> ko, in contrast to the case of weak in- 
compressible MHD turbulence. Q At = 45°, F k is £~ 7 / 2 
and Ak drops off more steeply than fc~ 7 / 2 . At = 7.1°, Fk 
falls off more rapidly than k~ n l 2 and AAF and AFF inter- 
actions keep Ak ~ F#. The cascade of fast-wave energy to 
high frequencies was found previously by 1 2 1 ] in direct nu- 
merical simulations with strongly turbulent Alfven waves and 
slow waves. In contrast to this Letter, these authors found an 
isotropic k^ 1 ! 2 fast-wave spectrum. 

The phenomenology described above can be applied more 
generally. For example, if the z-component of the phase veloc- 
ity, v p h lZ , were initially positive for all the excited waves, and 



4 




k/k 



FIG. 2: Top panel: power spectra as a function of k at 6 = 45°. 
Bottom panel: power spectra as a function of k at 6 = 7.1°. 



if there were no mechanism for generating Alfven-waves with 
k z = and v p h, z < 0, then there would be no AAA interactions. 
In this case, FFF interactions would still transfer fast-wave 
energy to high frequencies, and AAF and AFF interactions 
would still cause Ak and Fk to become approximately equal at 
small 8, but the Alfven-wave energy would not be swept out to 
large k± by AAA interactions. For waves with v p h Z > 0, one 
would thus expect Fk to obtain a constant-energy-flux k~ 1/!2 
scaling for all 8 with Ak ~ Fk at small 6. As a second exam- 
ple, if the initial excitation were primarily in Alfven waves, 
as may be the case in the corona |3], and if Ak were quasi- 
isotropic at k ~ ko, then AAF and AFF interactions would gen- 



erate significant fast-wave energy at k ~ ko, and FFF interac- 
tions would subsequently transfer fast-wave energy to higher 
frequencies. As a final example, if 8v nns <C va but the Alfven 
waves at small \k z \ became strongly turbulent at k± larger than 
some transition wave number k tr , as in 1 12], then collisions 
between oppositely directed Alfven wave packets would still 
transfer Alfven wave energy at any k z to larger k±, but the 
cascade time for this process at k± > k tY would change from 
1a to a new value, T^, str « £ ± 2 ^ 3 The Alfven waves in 

most of A:-space and the fast waves would still be weakly tur- 
bulent because the linear periods of these waves would still 
be much shorter than the nonlinear time scales, and much 
of the weak-turbulence picture would still apply. In partic- 
ular, Fk would be ^ kr 1 ! 1 for > 45°, and Ak and Fk would 
be steeper than AT 7 / 2 at small 8 and large k since Xa, Iamv, 
and Iaf are <C If in that part of £-space. 
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